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Abstract 

The excess of high energy neutrinos observed by the IceCube collaboration might originate from baryon number 
violating decays of heavy shadow baryons from dark mirror sector which produce shadow neutrinos. These sterile 
neutrino species then oscillate into ordinary neutrinos transferring to them specific features of their spectrum. In 
particular, this scenario can explain the end of the spectrum above 2 PeV or so and the presence of the energy gap 
between 400 TeV and 1 PeV. 


Recently the IceCube Collaboration published the 
data on high-energy neutrinos collected between 2010 
and 2013, containing 35 candidate events in the energy 
range from 30 TeV to 2 PeV, which show an evident 
excess over the expected background of the events with 
£>60-100 TeV or so fl] . On the other hand, no events 
were observed in the gap between 400 TeV and 1 PeV 
while three most energetic shower events emerged at 
the end of the spectrum with energies between 1 -2 PeV 
where the atmospheric background is practically van¬ 
ishing. The spectrum is apparently cut off at energies 
larger than about 2 PeV. The gap in the energy spectrum 
is difficult to explain in known models of high-energy 
neutrinos of astrophysical origin. 

Here we present a model ||2l that may explain such 
a spectrum. It is based on the idea that dark matter 
of the universe emerges from a parallel gauge sector, 
with particles and interactions sharing many similari¬ 
ties with ordinary particle sector. Such a shadow sector 
would contain particles like quarks which form compos¬ 
ite baryons, as well as leptons and neutrinos which are 
all sterile for ordinary gauge interactions. Particularly 
interesting example is represented by so-called mirror 
world a, which has the particle and interaction con¬ 
tent exactly identical to that of ordinary sector, with the 
same gauge and Yukawa coupling constants. 

Taking into consideration also attractive possibilities 
for physics beyond the Standard Model related to super¬ 


symmetric (SUSY) grand unified theory (GUT), one can 
consider that at higher energies our physics is presented 
by SUSY GUT, e.g. S U{5) or S U(6) which breaks 
down to the Standard Model S U(3)xS U(2)xU(l) at the 
scale Me - 2 x lO'® GeV. Supersymmetry breaking at 
Msb ~ 1 TeV triggers the electroweak symmetry break¬ 
ing and the Higgs field gets the vacuum expectation 
value (VEV) v = 174 GeV. In this view, we assume that 
at higher energies also mirror sectors is presented by the 
identical SUSY GUTs, S t/(5)' or S t/(6)', which breaks 
down to its standard subgroup S U{3)' x S 1/(2)' x 1/(1)' 
at the same scale Me - 2 x lO'® GeV. However, fol¬ 
lowing refs. 0 a, we assume that the symmetry 
between two sectors is broken later so that the elec¬ 
troweak symmetry breaking scale v' in mirror sector is 
much larger than ordinary electroweak scale. Namely, 
if y' ~ 10^' GeV, the lightest shadow baryons have 
masses order few PeV, and they decay due to baryon 
violating GUT gauge bosons, with decay time compa¬ 
rable to the age of the Universe, producing energetic 
shadow neutrinos which then oscillate into active neutri¬ 
nos (with oscillation probablities ~ or so) transfer¬ 
ring their spectrum to the latterQlt is worth to note that 
the decaying dark matter model, with a fraction of dark 
matter of about 10 per cent decaying before of present 
epoch could reconcile the Planck collaboration results 


*For other type of decaying dark matter model see e.g. Ref. (6). 
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on the CMB measurements with low redshift astrophys- 
ical measurements 0. 

In other words, we consider a supersymmetric grand 
unification theory SU(5) x S 11(5)' or SU(6) x S 11(6)'. 
As discussed in ref. 0, our proposal can be more nicely 
realized in SUSY S U(6) theory Q which gives natural 
solution to the so called hierarchy and doublet-triplet 
splitting problems via the Goldstone boson mechanism, 
relating the electroweak symmetry breaking scale to 
the supersymmetry breaking scale Msb, and in addi¬ 
tion naturally explains the fermion mass spectrum. In 
both sectors the GUT symmetries are broken at the 
scale Me - 2 x 10'® GeV. Below this scale our sector 
is represented by the minimal SUSY Standard Model 
(MSSM) S U(3) X S U(2) x U(l) with chiral superfields 
of quarks q, - (u,d)i, and leptons 1, = (v,e),, 

g; (i - 1,2,3 is family index) and two Higgs super¬ 
fields h and h, described by the Yukawa superpotential 
W - + Yf-gid^h + Supersymmetry is 

then broken at the scale Msb ~ 1 TeV inducing also the 
electroweak symmetry breaking. 

As for parallel mirror sector, below the scale Mq ^ 
2 X 10'® GeV we have supersymmetric 81/(3)' x 
SU(2)' X U(l)' theory with the similar particle con¬ 
tent, quarks q'. - (u',d')i, leptons /' = (v',e')i^ 

e':', and two Higgs superfields h' and Ji'. At the scale 
Me mirror gauge coupling constants j i ^re equal to 
the ordinary gauge constants g 3 , 2 ,i, and coupling con¬ 
stants in the Yukawa superpotential W' - Y/l'.eJh' + 
Y^.q'.d'^'h' + Yf.q'.u"JY have exactly the same pattern as 
in W. We assume, however, that in mirror sector super- 
symmetry and electroweak symmetry are both broken at 
the scale of about 10" GeV. 

Hence, below GUT scale, the gauge coupling con¬ 
stants ^ 3 , 2,1 and g '^2 1 evolve down in energies in both 
sectors in the same way up to scales of about lO" GeV 
where the supersymmetry is spontaneously broken in 



Logjg (/i/GeV) 

Figure 1: RG running of gauge couplings below the GUT scale in 
ordinary and shadow sectors, ai = 14;: and a'- = g'^ I An. 


mirror sector (see Fig. [^l. However, below this scale 
ordinary sector still remains supersymmetric and con¬ 
stants evolve down by the renormalization group (RG) 
equations as in the MSSM, down to scale Msb ~ 1 TeV 
where supersymmetry is effectively broken. After that 
the Higgses h and h are not protected anymore by 
the supersymmetry and they get VEVs vi = vcosyS, 
V 2 = vsinjS, V = 174 GeV, which induce the elec¬ 
troweak symmetry breaking and generate the fermion 
masses. The masses of lightest fermions, = 0.5 MeV, 
niu - 3 MeV and mg ^ 5 MeV respectively for electron, 
up-quark and down quark, are related to smallness of the 
Yukawa constants of the first generation. At the QCD 
scale A ^ 200 MeV gauge interactions of S U(3) be¬ 
come strong and confine the quarks into baryons, with 
lightest ones being proton and neutron with masses of 
about 1 GeV and spin 1/2. 

Let us consider now parallel mirror sector where su¬ 
persymmetry is broken at the scale Mgg ~ lO" GeV, 
supposedly due to non-zero F or D terms of some aux¬ 
iliary fields. Hence shadow scalars, including squarks 
and sleptons as well as the gauginos and the Higgs dou¬ 
blets h' and h' all acquire soft masses order Msb- Re¬ 
spectively shadow Higgses can get VEVs Vj = v' cos/S' 
and Vj = v' sin/S', which break mirror electroweak sym¬ 
metry at the scale v' < Mgg|^ Therefore, the masses 
of shadow fermions are rescaled, modulo renormaliza¬ 
tion factors order 1, by a factor ( - v'/v with respect to 
ordinary fermion masses. Namely, taking v'/v - 10^ 
and assuming tan/?' = tan^S, by the RG running of 
gauge and Yukawa constants from the GUT scale down 
in energies, one obtains for the masses of lightest mirror 
fermions Me - 0.4 PeV, Mo ^1.1 PeV and My ^1.9 
PeV, where capital letters E, D, U denote respectively 
the shadow electron e' , down quark d' and up quark u'. 
(Notice, that in mirror sector up-quark becomes heavier 
that in down-quark due to the difference in the RG run¬ 
ning of the Yukawa constants 0 .) One the other hand, 
the RG evolution of the S U(3)' gauge constant g'^ shows 
that the shadow QCD scale becomes A' ~ 100 TeV 
(c.f. A ^ 200 MeV in ordinary QCD), as it is shown 
on Pig. [T] Therefore, Mj/,M£) » A' and the shadow 
QCD looks like a rescaled version of our QCD but with¬ 
out light quarks, containing only the heavy quarks like 
c and b. In fact, Mj/ and Mo are larger than A' by about 
the same factor as the ordinary beauty and charm quark 
masses, and 111 ^, are larger with respect to ordinary 
QCD scale A. 


^Interestingly, if supersymmetry breaking is transferred to our sec- 
tor via gravity or other Planck scale mediators, this would nicely ex¬ 
plain ordinaiy soft masses order Msb ~ ^SB ~ ^ 
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As far as in shadow sector up quark U is heavier than 
down quark D, the lightest shadow baryon should be 
shadow baryon of spin 3/2, consisting of three down 
quarks D and having mass ~ 3Mo = 3.3 PeV. All 
states, containing up quark U, will be unstable against 
weak decays, U DEV. As for mesons, the light¬ 
est pseudoscalar is shadow neutral pion if consisting 
of DD, with mass Mq x 2Mo - 2.2 PeV, while the 
lightest vector meson p^iDD) is slightly heavier than n^. 
(Another neutral pion consisting of UU becomes much 
heavier, with mass of about 3.8 PeV.) Charged Pion 
as well as p*-meson consisting of DU will have mass 
M ^ Mu + Mo - 3 PeV, with a bit heavier than . 
All pseudoscalar and vector mesons have excited states 
with mass gap order A' between the levels, just like cc 
or bb states in our QCD. 

Now we come to the role of baryon violation and pro¬ 
ton decay which is fundamental prediction of the GUTs. 
The heavy gauge bosons of S U(5) with baryon violat¬ 
ing couplings between quarks and leptons induce the de¬ 
cay of the lightest ordinary baryons (proton, or neutron 
bound in nuclei), with lifetime Tp ~ ~ 

10^' Gyr or so, where ac is gauge coupling constant at 
the GUT scale Mq ~2x 10'® GeV El. In the shadow 
sector, the similar couplings of GUT gauge bosons 
should destabilize the shadow A baryon. However, tak¬ 
ing into account that the latter is much heavier than the 
ordinary proton, Mf^lnip ~ 10 ®, its lifetime must be 
about 30 orders of magnitude smaller than the proton 
lifetime. Hence we get ta ~ ~ 100 Gyr or 

so, comparable to the age of the Universe to - 14 Gyr. 

The principal decay mode of A baryon is in vec¬ 
tor mesons, A“ Pa + where generically is 

a superposition of shadow neutrino flavor eigenstates 
V^^^. Each decay produces monoenergetic neutrinos 
y/. with energies £; = - MflM^), where Mq 

is the mass of p“ meson and Mi, 2 ,... are the masses of 
its excitations. In Fig. 2(a) the spectrum of neutrinos 
produced by decay of galactic dark matter is shown by 
sharp peaks (solid blue) for Mq, Mi , M2 respectively be¬ 
ing 3.0,3.1 and 3.2 PeV. Due to close degeneracy be¬ 
tween the masses of A-baryon and p^ mesons, the neu¬ 
trino energies are £, <si 1 PeV while the most of initial 
energy = Ma is taken away by vector mesons p7. 

However, vector mesons readily decay into mirror 
pion and photon, p“ ^ tt + 7 ', and subsequent decay 
of the pion produces the neutrino once again (solid red 
curves in Fig. 2(a)). Shadow has two decay modes, 
two body 7 r“ —> e'v' and three body —» ifie'V^, where 
TT®, j - 0 ,1,... are the basic shadow pion and its excited 
states. Interestingly, the 2-body and 3-body branching 


Mirror neutrino fluxes 
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Figure 2: (a) Shadow neutrino spectrum produced by A-barion decay 
and subsequent prompt decays of shadow pions (b) Spectrum of VHE 
neutrino events at IceCube as predicted in our model. Magenta line 
shows contribution of the neutral cuiTent interaction and its integral 
corresponds to 1 event in the energy interval 100 TeV — 2 PeV. 


ratios are comparable. This fact is intimately related to 
the value A' ~ 100 TeV Two body decay produces 
neutrinos with a narrow energy spectrum concentrated 
around Ma/2 ^1.6 PeV, while the three body decay, 
due to smaller phase space, produces less energetic neu¬ 
trinos with a wide spectrum extending up to the value 
£max ^M_-Mo-Me^Mo-Mo-Me^ 0.4 PeV. 

Fig. 2(a) shows the final spectrum of shadow neu¬ 
trinos and v', including the neutrinos produced by 
the decay of dark matter in the galactic halo, and ex- 
tragalactic neutrinos produced by the decay of cosmo¬ 
logical dark matter at large redshifts. For definiteness, 
in Fig. 2(a) we assume that dark matter entirely con¬ 
sists of shadow baryons A, i.e. the fraction /a = 1, 
and take the decay time ta as 10 times the Universe 
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age, ta = The fraction of extragalactic neutri¬ 

nos strongly depends on the decay time ta, it increases 
with Ta decreasing. In correspondence, the energy gap 
becomes less pronounced and at ta < Q.lt^ or so it 
practically disappears since cosmological contribution 
from high redshifts become dominant. In this case dark 
matter cannot be entirely from mirror sector, the later 
can constitute only a smaller fraction and other type of 
stable dark matter should be also invoked^ 

The shadow neutrinos may have mixing with ordi¬ 
nary ones m. We assume that operators mixing 
and states respect a conservation of a combined 
lepton number L - L- L' d, and that all operators for 
neutrino masses are suppressed by the Planck scale Mpi 
Q. Hence, for neutrino masses and mixing we consider 
the following operators (family indices are suppressed): 


—^llhh 5 —H- ll'hh' 

Ml^ Ml, Mpi 


( 1 ) 


where x and x are gauge singlet chiral superfields with 
lepton numbers L - —2 and 2 respectively, with VEVs 
(t) - (t) - P breaking the lepton number. Alterna¬ 
tively, X and X can be promoted as flavon sextet and 
anto-sextet fields of the flavor symmetry SU(3)h im 
between three families assuming that it is a common 
gauge symmetry between two sectors IfT^ . This would 
give an certain guideline for obtaining predictive pat¬ 
terns for active and sterile neutrino masses and mixing. 

Here the first two terms give the Majorana masses re¬ 
spectively to ordinary neutrinos ^ and their shadow 
(sterile) partners v'^ ^ while third term induces the mix¬ 
ing (Dirac) terms between active and shadow neutrinos. 
Then total 6x6 mass matrix of v and v' states reads 113 : 


My 


niy niyy! 

mly, niy 


AT \ 

M \ AT^2 j 


( 2 ) 


where (, - v' jv and T = fi/Mpp Therefore, taking e.g. 
T ~ 10“^ and constants A,D ~ 10“^, we see that shadow 
neutrinos acquire masses ~ 10 keV and active sterile 
mixing angles are ~ 10“"^. 

Fig. 2(b) shows how the spectrum of shadow neutri¬ 
nos shown on Fig. 2(a) and transferred to ordinary neu¬ 
trinos via active-sterile mixings will be seen by the Ice- 
Cube. Here the effective areas for the neutrino detection 
by IceCube and characteristic error bars in estima¬ 
tion of neutrino energies (of about 13%) are taken into 


account. The obtained spectrum of events indeed look 
very much like the spectrum observed by the IceCube 
m. In addition, this model can has specific predictions 
for the flavor content of these high energy neutrinos in 
different parts of the spectrum lfT4l . 

The validity of our model can be tested with increas¬ 
ing statistics by the IceCube collaboration. 
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